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Abstract 

The article presents some aspects on the use of computer in teaching general rela- 
tivity for undergraduate students with some experience in computer manipulation. 
The article presents some simple algebraic programming (in REDUCE+EXCALC 
package) procedures for obtaining and the study of some exact solutions of the 
Einstein equations in order to convince a dedicated student in general relativity 
about the utility of a computer algebra system. 

1 Introduction 

Teaching general relativity is a very difficult task not only for the "teacher" 
but also for the students. But why ? 

General relativity (GR) is not only a theory of gravity; it is a theory of 
the structure of space and time, and hence a theory of the dynamics of the 
universe in its entirety. Thus, the theory is a vast edifice of pure geometry, 
indisputably elegant, but of a great mathematical difficulty especially for 
undergraduate students (among others...). But why undergraduate students? 
Because, in our department, (considering GR a necessity for the medium 
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graduate physics people) we make efforts to introduce a course of GR at the 
Ill-rd year level. After some years of experience in teaching GR at graduate 
students, specialized in theoretical physics, we can say that after some weeks 
of introducing the must important tools of differential geometry, starting 
with the physical problems of GR, an important part of our students were 
almost discouraged. At an undergraduate level, with unspecialized students, 
this risk is much greater. 

Thus, algebraic programming systems (like REDUCE - [|J) which con- 
tain differential geometry packages can become a very important tool for 
surpassing these difficulties. With the computer, the student can learn very 
fast, and in an attractive manner, the important notions of differential ge- 
ometry, tensor calculus and, of course, the exterior calculus (with EXCALC 
for REDUCE, for example - @). As an example, he can make, after some 
simple computer manipulation, the long and unattractive (for the "man- 
ual" student) calculations from the Riemannian geometry, with metrics and 
Christtoffel symbols, etc. Thus, the first section of this article illustrates how 
we can use the EXCALC package in teaching Riemannian geometry. 

But having such a powerful tool for calculus in differential geometry, we 
can try to introduce some aspects of GR using computer algebra systems. For 
example, the Schwarzschild solution can be easy obtained on the computer, 
and also we can generate the Reissner-Nordstrom solution as exact solutions 
of the Einstein equations. These aspects are presented in the second section 
of the article. 

The last section of the article presents the use of computer algebra in 
finding and treating other exact homogeneous solutions of the Einstein equa- 
tions with cosmological constant (de Sitter and anti de Sitter metrics more 
precisely) . 

As a conclusion, we consider the use of computer as an important tool for 
teaching general relativity. During the last two years we have experienced 
several packages of procedures, (in REDUCE + EXCALC for algebraic pro- 
gramming and in Mathematica or Maple for graphic visualizations) which 
fulfill this purpose. Even when the students were real beginners in computer 
manipulation we have obtained visible good results, in approaching several 
topics of differential geometry and of course, in general relativity. 
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2 Differential geometry in EXCALC 



The program EXCALC (0) is completely embedded in REDUCE, thus all 
features and facilities of REDUCE are available in a calculation. 

EXCALC is designed for easy use by all who are familiar (or want to 
became) with the calculus of Modern Differential Geometry. The program is 
currently able to handle scalar-valued exterior forms, vectors and operations 
between them, as well as non-scalar valued forms (indexed forms). 

Geometrical objects like exterior forms or vectors are introduced to the 
system by declaration commands; therefore zero-forms (functions) must also 
be declared. Also, specific operations with geometric objects are available in 
EXCALC like : exterior multiplication between exterior forms (carried 
out with the nary infix operator A (wedge)), partial differentiation (is 
denoted by the operator @), exterior differentiation of exterior forms 
(carried out by the operator d), the inner product between a vector and an 
exterior form (represented by the diphthong _ | (underscore or-bar)), the Lie 
derivative can be taken between a vector and an exterior form or between 
two vectors (represented by the infix operator | _), the Hodge-* duality 
operator (maps an exterior form of degree K to an exterior form of degree 
N-K, where N is the dimension of the space). It is possible to declare an 
indexed quantity completely antisymmetric or completely symmetric. Some 
examples : 

PFORM U=1,V=1,W=K; ^declaration of some forms 

(3*U-A*W)~(W+5*V)~U; 

A*(5*U"V"W - U"W"W) 

@(SIN X,X); % partial differentiation 

COS(X) 



PFORM X=0,Y=K,Z=M; 

D(X * Y) ; °/ exterior differentiation of a 

X*d Y + d X~Y % product of two forms 

D(X*Y~Z) ; 
K 

( - 1) *X*Y"d Z + X*d Y~Z + d X~Y~Z 



PFORM X=0,Y=K,Z=M; TVECTOR U,V; 
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U_ I (X*Y~Z) ; 7. inner product 

K 

X*(( - 1) *Y"U_|Z + U_|Y~Z) 



PFORM Z=K; TVECTOR U; 
U |_ Z; 



7o Lie derivative 



U_|d Z + d(U_|Z) 

A metric structure is defined in EXCALC by specifying a set of basis 
one- forms (the coframe) together with the metric. The clause WITH MET- 
RIC can be omitted if the metric is Euclidean and the shorthand WITH 
SIGNATURE <diagonal elements> can be used in the case of a pseudo- 
Euclidean metric. The splitting of a metric structure in its metric tensor co- 
efficients and basis one-forms is completely arbitrary including the extremes 
of an orthonormal frame and a coordinate frame. Examples (0) : 

COFRAME 0(T)=D T, X=D X 

WITH SIGNATURE -1,1; %A Lorentz coframe; 

COFRAME E R=D R, E PH=D PH %Polar coordinate 

WITH METRIC G=E R*E R+R**2*E PH*E PH; 7„basis; 

The frame, dual to the coframe defined by the COFRAME command 
can be introduced by FRAME <identifier>. This command causes the dual 
property to be recognized, and the tangent vectors of the coordinate functions 
are replaced by the frame basis vectors. 

The command RIEMANNCONX is provided for calculating the connec- 
tion 1 forms. Example : calculate the connection 1-form and curvature 
2-form on S(2) (displaying only the nonzero results) : 

COFRAME E TH=R*D TH,E PH=R*SIM(TH) *D PH; 
RIEMANNCONX OM; 

0M(K,-L); 7„Display the connection forms; 



PH 



PH 

:= (E *COS(TH))/(SIN(TH)*R) 



NS 



TH 
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TH PH 
NS := ( - E *COS(TH))/(SIN(TH)*R) 

PH 



PFORM CURV(K,L)=2; 

CURV(K,-L) :=D OM(K,-L) + OM(K,-M) "OM(M-L) ; 

%The curvature forms 



PH 



TH PH 2 
( - E ~E )/R %it was a sphere with 

°/ radius R. 



CURV 



TH 



TH 



TH PH 

(E ~E )/R 
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CURV 



PH 



3 General relativity on the computer. 
Schwarzschild solution 

The students in our Faculty of Physics are, generally speaking, well trained in 
practical computer manipulations. There is no semester without at least one 
course with labs in the computer room. But when we invited our students 
to come in the computer room, to learn something about general relativity 
with the computer, it was a general surprise, because they considered (until 
now) the computer as a tool for hard numerical computations. They do 
not know almost anything about computer algebra - [[|. 

It is not necessary to use sophisticated procedures, with large and com- 
plicated metric statements (which are almost impossible to calculate by hand 
in a civilized time of teaching) in order to convince a dedicated student in 
general relativity about the utility of a computer algebra system. It is enough 
to run a simple program like (||) : 

pform psi=0; f domain psi=psi(r); 



cof rame 



o(t) 
o(r) 



= psi 

= (1/psi) 



* d t, 

* d r, 



% Schwarzschild 
% metric 
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o(theta) = r * d theta, 

o(phi) = r * sin (theta)* d phi 
with signature 1,-1,-1,-1; frame e; 

to introduce a Schwarzschild type metric in spherical coordinates (r, 9, ip). 
This means that in classical notation the interval is 

ds 2 = V 2 dt 2 - -^dr 2 - r 2 (d9 2 + sm 2 6 dip 2 ) (1) 

^ being a function \l/ = \l/(r) and differ from a Minkowski one by the a new 

"unknown" function = \J\ + unknown(r) which must be determined 
from Einstein equations : 

pform unknown=0; f domain unknown=unknown(r) ; 
psi := sqrt(l + unknown); 

Now comes the most important part of the procedure : the calculation 
of the components of Einstein tensor (einstein3) via the Riemann or Levi- 
Civita connection 1-form r y ' - christl and the curvature 2-form R % i - curvl 

pform chrisl (a,b)=l , curv2(a,b)=2, einstein3(a)=3; 
antisymmetric chrisl, curv2; 

riemannconx christl; chrisl(a,b) := christl (b, a) ; 

curv2(a,b) := d chrisl(a,b) + chrisl(-c,b) ~ chrisl(a,c); 

einstein3(-a) := (1/2) * curv2(b,c) ~ #(o(-a) ~o(-b) ~o(-c)) ; 

The last of the above program lines just defines the Einstein 3-form which 
appears in the Einstein equations. Those who prefer the coordinate compo- 
nents form (thus the Einstein tensor Gij = Rij — \gijR) can use the next line 
to "pick-up" these components as : 

pform Ein(i,j)=0; 

Ein(-i,-j) :=e(-i)_|#einstein3(-j) ; 

A typical component (G^) of the output of einstein3 reads : 
T R THETA 

- (0 "0 ~0 *(@ UNKN0WN*R + 2*@ UNKNOWN) )/ (2*R) 

R R R 
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or : 




1 d 2 unknown 1 dunknown 




) 



& A o T A of 



Requiring the coefficients to vanish yields a second order differential equa- 
tion for the function unknown. Trying unknown = a * m/r * *n, after 
using SOLVE package (i.e. the EXCALC command solve(einstein3(- 
phi), unknown); ), we obtain n = —1 and : 

unknown := - alpha * m/r; 

where "m" is the mass and alpha a constant coefficient to be determined by 
physical considerations (link to the Newtonian theory, for example). 
Finally, evaluating the ip function (psi := psi), we obtain : 



which is the typical form of the Schwarzschild metric (@|-0) identifying 
a = 2 by physical considerations. 

From now one, it is possible to study, in a similar way more complex 
situations, like Reissner-Nordstrom metric (starting, of course with the above 
Schwarzschild one - this example is presented in detail in 0). The teacher 
can select more exact solutions of Einstein equations in order to complete 
the education of his students. Also, algebraic programming can be used to 
present (in a very fast manner) the canonical version of general relativity 
(|0]) or the post-Newtonian approximation ([||). 

The next section is dedicated to more examples which can be used in the 
teaching process of general relativity. 

4 Other examples 

According to the literature (see, for example |3J or ||) two of the homoge- 
neous solutions of the Einstein equations with cosmological constant are the 




or, in ([I]) we have 
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de Sitter and anti-de Sitter metrics, having the line element written as 



ds 2 = e n9 dt 2 - e kf dx 2 - e n9+kf (dy 2 + dz 2 ) (3) 

where / and g are two unknown functions of x and t variables which we shall 
determine by imposing that the vacuum Einstein equations with cosmological 
constant to be fulfilled : 

Gy = Rij - -gijR - Agij = (4) 

Also in the above equation (|3|) n and k are two constants introduced in order 
to control the two solutions : for n = and k = 1 we shall obtain the de 
Sitter solution and for n = 1 and k = 1 we have the anti-de Sitter solution. 
Thus we can use a program sequence for introducing this metric as : 

pform f=0,g=0; fdomain f =f (x,t) ,g=g(x,t) ; 
cof rame 

o(t) = exp(n*g/2) * d t, 

o(r) = exp(k*f/2) * d x, 

o(theta) = exp((n*g+k*f )/2) * d y, 
o(phi) = exp((n*g+k*f )/2) * d z 
with signature 1,-1,-1,-1; frame e; 

Now we must run the sequences which calculates the Einstein tensor 
components with cosmological, namely : 

riemannconx chris; 
pform chrisl(a,b)=l; 
pform curv2(a,b)=2; 
pform einstein3(a)=3; 
antisymmetric chrisl; 
antisymmetric curv2; 
chrisl(a,b) :=chris(b,a) ; 

curv2(a,b):= d chrisl (a, b)+chrisl (-c ,b) "chrisl (a, c) ; 
einstein3(-a) :=(l/2)*curv2(b,c)" #(o(-b)"o(-c)~o(-a))-Lam*#o(-a) ; 
pform Ein(i,j)=0; 

Ein(-i ,-j ) :=e(-i) _ | (#einstein3(-j ) ) ; 
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where Lam represents the cosmologic constant A and Ein(-i,-j) represents 
the components of the Einstein tensor Gij with cosmological constant in 
coordinate frame. A typical output is, for example : 

F*K + G*N F*K 2 2 

EIN := ( - 4*E *LAM + 3*E *@ F *K 

T T T 

F*K F*K 2 2 

+ 4*£ *@ F*@ G*K*N + E *@ G *N 
T T T 

G*N G*N 2 2 

- 4*£ *@ F*K - E *@ F *K 

XX X 

G*N G*N 

- 4*£ *@ F*@ G*K*N - 4*E *@ G*N 

XX XX 

G*N 2 2 F*K + G*N 

- 3*E *Q G *N )/(4*E ) 

X 



which means, in "normal" form : 

2 



G» = e 



-ng 



C [dt +kn dtdt + A n U 



dx 2 4 \ <9x / dx dx dx 2 4 \dx 



-A 



Now assigning to the parameters the values n = and k = 1 the above 
component of the Einstein tensor) and the corresponding Einstein equation) 
will be : 



r 3 ( d n f 
Gtt = -Am ~ e 



dx 2 4 1 dx 



- A = 
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By inspecting this equation (0) we can observe that his general solution can 
be written as f(x,t) = fl(t) + f2(x) where 

V dt J 3 3 

d 2 fm _ 1 f2fa) 1 (df2{x)V 
dx 2 4 1 4 ^ da; / 

where C\ is a constant. Thus we can choose f(x,t) = fl(t) such that the 
above equation becomes : 

pform f l=0;fdomain f l=f 1 (t) ; f : =f 1 ; 
Ein(-t ,-t) ; 
2 

3*@ Fl - 4*LAM 
T 



4 

as computer output and 

This equation has an obvious solution of the form : 

fl = f=^V3At + C 1 

which fulfill the Einstein equations with cosmological constant (all the com- 
ponents os Ein(-i,-j) 0-form are zero). Then the line element becomes : 

ds 2 = d t 2 - e lV3At +Cl ^2 + dy 2 + ^ 

which is, of course the well-known De Sitter solution. For the anti-De Sitter 
solution we must choose the parameters as n = 1 and k = 0. Thus the 
interesting Einstein tensor components become : 
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n:=l;k:=0; 
Ein(-t ,-t) ; 



G G 2 G 2 

- 4*£ *@ G - 3*E *@ G - 4*E *LAM + @ G 
XX X T 



4*E 



Ein(-x,-x) ; 



G 2 G 2 

3*E *@ G + 4*E *LAM - 4*0 G - @ G 
X T T T 



4*E 



as computer output, or in "normal" transcript : 



G, 



G tt = : 

3 /c^ 
r ~ 4 lar 



i v 5t . 

2 



cte 2 4 ycto^ 

d 2 g \ (dg\ 
W + 4 [di 



+ A = 



Solving, in a similar way as in the De Sitter case, the second of the above 
Einstein equations, we have to impose g(x,t) := gl(x). Thus we obtain a 
solution as : 



g = gl = -V^3Ax + C 2 
where C 2 is a new constant. We have the line element as 

ds l = _ dx 2 + efv^H-Cx ( rff 2 _ dy 2 _ d ^ 

which is the anti-De Sitter solution. 
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The REDUCE+EXCALC procedures we presented here are also suitable 
for a series of other applications during the relativity lab classroom (...). For 
example the teacher can propose to the students to check if other metrics 
fulfill the Einstein equations. They can use the above program sequences with 
minor modifications. As an exercise we propose here the Ozsvath metric (||) 

ds 2 = dx 2 + e 2 ^T^> x (dy 2 + 2dudv) + f e -^ Z] ^ x (-2y/2dy + fe^V^^dvj 

which is also an homogeneous solution of vacuum Einstein equations with 
cosmological term. Here a special care must be taken to the changing of the 
signature of the metric and to the nondiagonal terms of the metric. 
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